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Abstract 

A more rigorous treatment of the Schwarzschild metric by making 
use of the energy-momentum tensor of a single point particle as source 
term shows that 

5oo = -{1 - ^ - *S^(0(r) - l)}exp[2(0(r) - 1)] 

9rr = {1 - - *g£(*(r) - I)}" 1 

The existence of a discontinuity at r = leads to an infinite repulsive 
force that will change the ultimate fate of a free fall test particle to a 
bouncing state. 



PACS: 04.20.Jb , 04.70.Bw 
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We intend to show that at r = there is a mathematical inconsistency in 
derivation of the Schwarzschild metric which is the static solution of Einstein 
field equations for vacuum space around a point mass M[l] . It is a simple 
case that is propounded in any related text in GR. Usually similar approaches 
to the problem are followed [2-6]. We use the notations and derivations of 
Weinberg [2] for the sake of defmiteness. 
The standard form for general isotropic metric is : 

dr 2 = B(r)dt 2 - A(r)dr 2 - r 2 d6 2 - r 2 sin 2 6d V 2 (1) 

The components of Ricci tensor for the metric (1) are: 

B" B'X B\ A' 
Rrr -2B~4B { A + W-VA (2) 



R w = sin 2 6 Ree (4) 



B" B'XB'B' 
R - = -2A + U^A + ^-VA (5) 



R^u = for n^v 

Prime means differentiation with respect to r. The applied field equation for 
empty space is 

Rfw = o (6) 

Using (6) we may conclude that ^ + = , and by inserting (2) and (5) 
in this relation we get 



2 



1 .A' b\ 

Integration of (7) with respect to r and imposing Minkowski metric as r 
tends to infinity for boundary condition gives 

A=i (8) 
Now by arranging (2) ,(3), (6) and (8) we obtain 

Ree = -1 + rB' + B = (9) 



B B 

R rr = — + — = (10) 
IB rB K J 

Integrating (9) with respect to r gives 

rB = r + const. (11) 

For fixing the constant of integration in (11) we recall that at great distances 
from mass M,the component of gu = —B must approach to — 1 — , where 
$ is the Newtonian potential — Hence it is equal to — and our final 
result is 

B=l-*^ (12) 

c l r 

According to the correspondence principle general relativity must agree 
on the one hand with special relativity in the absence of gravitation and on 
the other hand with Newtonian theory of gravitation in the limit of weak 
gravitational fields and low velocities. The corresponding field equation of 
this problem in Newtonian gravity is the Poisson equation 

V 2 $ = AnGM5(7) 

It is remarkable to notice that this equation is defined in the whole space 
even at r = and its solution which is proportional to inverse of r satisfies 
the field everywhere even at r = 0. Thus we have right to expect that the 
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general relativity field equations which are supposed to be yielded to this 
equation at weak field limit must be well defined at r = [7,8]. 

Since the singularity at r = is not a hypothetical case and it is the 
inevitable end state of gravitational collapse of a massive body , the external 
gravitational field of our point particle must somehow be traced to the source 
located at r = 0. We consider this classical problem not in the belief that 
it is the final word , but rather a useful start which quantum phenomena 
will no doubt modify. It can then be used as a background for the more 
advanced theories. Such advanced theories are important , but they have 
achieved their greatest successes when the corresponding classical problem 
has been fully understood. 

The important point is that for this case does not vanish everywhere, 
it is nonzero at r exactly equal to zero. For finding the T^ u of a mass point 
M located at x 1 we may use the action [9] 

'» = L M-lrW ™ (13) 
is defined to be 

SI M = \J d 4 x^/gJx)T^(x)8g^(x) (14) 
A straight forward calculation gives 



1„, f +OD . r , , .sdxKp) dxi(p), i dxHp) dxUp) . , , 

(15) 

This is of the form (14) with 



T Xk (x) = g^(x)M dr^^8\x - Xl ) (16) 

if we take / d i xS 4 (x — xi) = 1 as is taken by Weinberg. It would be more 
concrete if we take / ^fgd A x5 4: (x — X\) — 1 which in this case we come to 

m\ K / \ a, f + °° r r / \ dx% dxi,_i dx\(p) dxi(p) r4 , , XN 
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For a partcile located at the origin we have x% = 0) , and in a frame 
at rest with respect to the particle the components of the four velocity are 

v} = , u* = 5-5 , Ut = -53 (18) 



Thus the (OO)-component of (17) is 

dp[B(0)]-^S 4 (x - xi) 
-00 dp 

= Mp^n^l) (19) 

To find out the form of 5 A (x — x\) it is more convenient to work in the 
coordinate system x^ = (x°,r,q = cos9,ip). In this system we have g qq = 
and y/g = r 2 V ' AB. Thus we may write 
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J r 2 VAB5 4 (x - Xl )d A x = 1 (20) 
(20) will be satisfied if we take 

where we have J +o ° 5{r)dr = |. 
Inserting (21) in (19) gives 

!«(*) = M ^ (r) (22) 
2irr 2 y/ABB(0) 

We should notice that T tt = c 2 T 00 and T rr = BT tt . Thus we have 

, Mc 2 5(r) . . 

T tt (x) = v ; (23) 



2irr\/ABB(0) 

and 



T rr = M*M (24) 
2vrr 2 V^(0) 
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Using (24) we may write 

J T TT 4nr 2 VAdr = Mc 2 (25) 
as we should have expected. Also we have 



Then T = g^T^ is 



Mc 2 5(r) rB 
2nr 2 \j AB(0) 1 ' 



Now the components of which are defined as S^ u = T^ v — \g^ u T are: 



Mc 2 8(r) 

Stt -^^ B iAB{0) ( ] 



_ Mc 2 5(r) VAB 



See = ^fi± (30) 



4 V 



AB(0) 



S w = See sin 2 (31) 



«V = fi ^ v (32) 
The rigorous treatment of the problem is to solve the following equations 

RfXV = ~^~Sfxv (33) 
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where the components of S^are given by (28)- (32). Using (33) and (2)- (5) 
we obtain 

2GM5{r) [IF 
rr ~ 2B ~ IBM + ~B> ~ VA ~ 1(0) ( ' 



r A' B' 1 2GM5(r) I B 
R °° = - 1 + 2A ~A + 13 ) + A = ^ABW) (35) 



R w = sin 2 9 R ee (36) 



~ ~2A 4~A ~A ~B ~ 7a~ 6V - *V AB(0) (3?) 

It can be checked easily that the usual solutions AB = 1 and B = 1 — 2 ^M- 
do not satisfy (34)- (37). This point has been recognized by Narlikar too 
[10]. He says when we try to solve for a point mass gravitational field the 
condition AB = 1 yields to an inconsistency at r = because this gives 
T° = T 1 ! while for a static problem we need T 1 ! = 0. He argues that it 
would be tempty to take the easy way out of the problem by taking that the 
Schwarzschild coordinates are in appropriate at r = 0. 

Using (34) and (37), + gives 



J_A B AGM5(r) B 

rA A B c 2 r 2 \j AB(0) 1 } 



and 



A B _ AGM8(r) / AB 

A B cV V B(0) 1 J 



Now we may combine (39) and (35) to write 



rA 2GM5(r) / B 1_ 2GM5(r) / B 

~A 2 + c 2 ^AB^j + A~ c 2 VAB(0) ( ^ 
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which gives 



Integrating (41) with respect to r from +00 to r and taking that 

= 77= ■ y' elds 



r 4GM(6(r) - 1) 
r + — + y ; ; = const. , (42) 



1 x > 
where 9(x) — \ \ x — . 

x < 

For finding the integration constant in (42) we may impose this fact that 
at large distances we have AB = 1. Then we get 

— r + rB + = const. , (43) 



and at this range B is equal to 1 — ^qr- By inserting this value for B in (43) 
the constant in (42) and (43) will be found to be equal to — ^3^. Putting 
this in (42) results that 

(44) 



2GM n _|_ 0(r)-l 

Taking the limit of r — > 0, (44) gives 



1_ 2GM[1 + 

c r \/^(0) 



r-0 2GM 

or 

/ c 2 r \ 



A ^) = ^4^4m (45) 



A(0) = lim — — (46) 



Considering (46), we may write (44) as 

1 

~-i0ir) - 1) 



1 2GM 8G 2 M 2 ^ i 7 (47) 
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because the correction term is merely effective at r = 0. 
Then we have 

A> 2GM _ 16G 2 M 2 (o( r \ _ -i \ i 8G 2 M 2 8(r) 

f±_ c 2 r 2 c 4 r3 \ a V ) c 2 r 2 (ao\ 

A 1 - m - *g£(e(r) - 1) 1 } 



Putting (48) in (39) and using this fact that at r.h.s of it \J -j^j indeed is 
A(0) and (46) yield 

B' m- + y%£-{9{r) - 1) 



(49) 



B 1 - ^ - *g£{9{r) - 1) 

Integrating (49) with respect to r from +oo to r and imposing B(+oo) = 1 
gives 

B = {l--± r - ^rV(r) - 1)} exp[2(^(r) - 1)] (50) 
Using (47) and (50) we have 

AB = exp[2(0(r) - 1)] (51) 



Finally our rigorous result for Schwarzschild metric is: 



ds 2 = c 2 {l-^-^^(^)-l)}exp[2(%)-l)]^ 

' dr 2 -r 2 (d9 2 + sin 2 6 dip 2 ) (52) 



2GM 8G 2 M 2 



(0(r) - 1) 



It can be easily verified that this metric satisfies the field equation exactly 
at each point even at r = 0. In a mathematical rigorous treatment one may 
truly expect that the theory of distribution should be considered [11-15]. 
This is not a possible task at this stage because distribution in curved space- 
time is not well known. Here theta and delta functions have been considered 
in the same manner which are commonly being used in physics literature. 
This work is a step toward a perfect treatment in this sense. Another point 
to notice is that in Newtonian gravity the potential is a continuous function 
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even at r = and the force is always attractive. In this result there exist a 
discontinuity in the metric at r = 0. The neighborhood of r = that is where 
r << 2M is a domain in which the field is very strong and therefore there is 
no weak field limit corresponding to it. This phenomenon is a pure general 
relativistic effect which has no analog in Newtonian gravity. The existence 
of this discontinuity means that we have an infinite repulsive force at r = 
(see appendix). For more clarification we will find the tidal forces and the 
scalar invariant v \ p R ll v>,p for this metric in the next section. This repulsive 
force causes to change the ultimate fate of a free fall test particle from a 
collapse into intrinsic singularity to a bouncing state. The prediction of this 
phenomenon itself may be considered as a justification for this calculations 
which was rather detailed and the presented treatment for determining the 
form of the line element. 



The tidal effect results in an elongation of the distribution in the direction of 
motion and a compression of the distribution in transverse directions. The 
same effect occurs in a body falling towards a spherical object in general 
relativity. We can gain some idea of this by considering the equation of 
geodesic deviation in the form [16] 



for the spacelike components of the orthogonal connecting vector rf connect- 
ing two neighbouring particles in free fall. 
Let the frame ei a be defined in Schwarzschild coordinate as 



1 Tidal Forces 




R a bcd e a a v h v c e p d ^ = 



(53) 



Dt 2 



e 
ei 



a 



a 



a 



= 5-3(1,0,0,0) 
= A-i(0, 1,0,0) 
= r-^O, 0,1,0) 



(54) 



a 



= (r sine)' 1 (0,0, 0,1) 



and denote the components rf by 



r? = (r ] \r ] \r?) = (r ] r ,rf,rf) 



(55) 
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Using (1) the nonvanishing components of the curvature tensor are 



pt 


b" 

2B 




TDT 

-n trt 


b" 

2A 


R l ete 


rA 
~ 2 A 2 




R T Brd 


= R t ete 


TDt 


= R t ete 


Sin 2 6 


or 

xt (pr<p 


- R* 


R e tet 


. b' 

2rA 




Rftyt 


= R d tet 


R 9 r6r 


A 

~ 2rA 




w 


= R 9 r 9r 




= R^Oif 


9 Sin 2 6 


Rfe^e 





By using (56) and (54) , (53) gives the following 

D 2 rf B" 



Dt 2 2AB 



r] r = 



D 2 rf B' e n 

H V = 

Dt 2 2rAB 1 

D 2 rf B' 



Dt 2 + 2rAB 
Using (47), (50) and(51) we have 

B' GM 8G 2 M 2 



rf = 



+ 



, 1N ±G 2 M 2 . <J(r) 



2AB c 2 r 2 c 4 r 3 v ~ v 7 ' c 4 r 2 ' i 

By (46) we have, = ^^5{r) , then (60) becomes 

B' GM 8G 2 M 2 



2AB c 2 r 2 + c 4 r 3 



-(9(r) - 1) 



Next (34) gives 
B" 



B ,A B\ A 2GM IB 



2AB AAB A B' rA 2 c 2 r V AB(0) 



11 



Inserting (61) , (39) and (46) in (62) yields 

B" 2GM \§G 2 M 2 lnl . . GM 

i^)-l) + -^6(r) (63) 



2AB c 2 r 3 c A r A v v ' ' c 2 r 2 

Finally putting (61) and (63) in (57)-(59) leads to 

D 2 rf r GM 8G 2 M 2 /n/ . iV1 

—^—77" = { — — - — - — (0(f) - l)}v 65 

Dt 2 1 c 2 r 3 c 4 r 4 y y ' n ' y ' 

D 2 rf f GM 8G 2 M 2 /n/ . _ w 

(64)-(66) exhibit the usual tidal force of the Schwarzschild metric for r ^ 0. 
At r = the correction terms are dominant and have opposite sign to the 
first terms. 



2 Scalar Invariant 

Since the metric components have no cross terms and furthermore the cur- 
vature tensors have some algebric properties, the Riemann tensor scalar in- 
variant take a simple form as 

R a bcd R a bcd = 2Y,(g a aR b ah aY (67) 

By using (56) we may write (67) as 

R\ cd Ra bcd = (|^) 2 + 4(^) 2 + 1(1 - \f (68) 
Inserting (63) and (47) in (68) gives 

TT R^-A\ 2GM 16G2M2 (0(r) ^ + GM Mr)\ 2 
R bcd R a -4{— — ^j—(V(r) ~ 1) + ~^o[r)\ 

4 2GM 16C 2 M 2 . 8G 2 M% 

4 2GM 1QG 2 M 2 , n . . 1X12 
+7^{^ + ^i7^(^)-l)} 2 (69) 
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(69) shows that R a bcd R a bcd behaves like 48 %f 2 for r ^ 0, the same result of 
the Schwaezschild metric. At r = in addition of this term there are some 
correction terms which all of them diverge. 

3 Remark 

According to singularity theorems a condition for the existence of the singu- 
larity is that R ab K a K b > for every non-spacelike vector K. Now using the 
four velocity (18) as a vector we obtain in this case 



R „ V - GM6{r) - 6{r) (70) 

RabU u —\IabW) - (70) 

which is negative for r = 0. This means this problem does not satisfy the 
conditions of singularity theorems and they are not applicable in this case. 
This is the reason why the test particle has the chance to escape this curvature 
singularity. 

It is noticable that from (46) we may infer that A and B^ 1 approach to 
positive zero as r goes to zero. This means at r = the common convention 
of t as time coordinate and r as space coordinate remain valid. This is in 
agreement with our initial assumption that the point particle is located at 
r = for ever. 

We may conclude our discussion by expressing that there is no trace of those 
conceptual problems which was raised by Narlikar in this form of the metric. 



4 Appendix 

The geodesic equations of this metric are the following, where p is a parameter 
describing the trajectory. 

rfV A' fdrV r fd9\ 2 rsm 2 9 ( d<p\ 2 f dtV 

df + ZA\dp A\dp^ A~ \~dp~ + 2A\d^> "° " h 



d 2 6 | 2d0dr sin 2 9 ,. r = () ( _, ?) 
dp 2 rdpdp ° 
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d 2 ip 2 dip dr ^ ^ ^ d6 
dp 2 r dp dp dp 



(73) 



dH B' dr dt _ Q 
dp 2 B dp dp 



(74) 



Since the field is isotropic, we may consider the orbit of the test particle 
to be confined to the equatorial plane, that is 9 = |. The equation (72) 
immediately is satisfied and we can forget about 9 as a dynamical variable. 
Then equation (73) gives: 



d 2 ip 2 dip dr 
-rrr-\ r-— = or 



d 
dp 



dp 2 ' r dp dp 
Integrating (75) with respect to p leads to: 



In 



' dip s 
dp, 



lnr 2 







dip 
dp 



J_ 



where J is a constant. Dividing (74) by ^ gives 



d 

dp 



In 



'dt_' 

dp. 



+ lnB 



or 



dt 
dp 



1 

B 



(75) 



(76) 



(77) 



The related constant of integration has been absorbed into the definition of 
p. Inserting (76) into (71) we obtain 



d 2 r A' ( dr s 
dp 2 2 A \dp j 



J 2 



B' 



r 3 A 2AB 2 



(78) 



or 



d 

dp 



A(r) 



dr 



J 2 



dp J r 2 B 
By integrating (79) with respect to p we get 



(79) 
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. , . / dr \ 2 J 2 1 „ . . 

M*J (80) 

where E is a constant of integration and if the test particle is going to be 
moveless at infinity we must have E — 1. Equation (80) for a geodesic with 
dip — i.e. J = becomes 

^©'-i- 1 (81) 



Inserting (81) in (78) with J = gives 



(fr 1 / A B \ A 



dp 2 2AE V ^ 5 / 2yl2 
and inserting (39) in (82) gives 



(82) 



dp 2 2drVA7 (? ry fABB(0) 
and at last using (47) and (50) in (83) yields the final conclusion 

d 2 r GM 8G 2 M\ n/ . 1X .4G 2 M 2 n . r/ . 

^ = - *s - — {0{r) - 1] + '— - 2eWr) (84) 

In the r.h.s of (84) the first term exhibits the common attractive force and 
the other terms represent point interactions which act at r = 0. Indeed at 
the origin the third is the dominant repulsive force. 
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